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Abstract 

The phase "jumps" for sohtons interacting on a vortex filament, 
observed in experiments, have been unaccounted for since more than 
twenty years. Using exphcit formulas describing the interaction of two 
sohtons on a thin vortex filament in the Localized Induction Approx- 
imation we show that an appropriate choice of the parameters of the 
sohtons leads to large phase shifts. This result does not depend on 
the axial fiow along the filament. 
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It is well known that the equation describing the motion of a single thin 
vortex filament in an ideal fluid in the localized induction approximation 
P is integrable in the sense of the soliton theory [2]) and admits soliton 
solutions (known also as kinks or Hasimoto vortices) It is interesting 
that the equation itself has been derived as early as 1906 by da Rios and 
then rediscovered few times, see However, the experiments carried out 
some time later have shown that although some soliton-like patterns can be 
observed and their qualitative behavior agrees with the theory (j3]) but there 
are some discrepancies concerning quantitative details j6j . Let us mention the 
ratio of the group velocity and the phase velocity, the dependence between 
phase velocity, angular velocity and group velocity, and last but not least: 
the very large "phase jump" during head-on colhsions of two solitons [71IH]. 

Then, the approximation has been improved by taking into account the 
axial flow (the flow inside the vortex core, along the filament) [H]. The 
kinematics of the vortex filament remained integrable but the equation is 
now more general (it contains one additional free parameter corresponding 
to the intensity and direction of the axial flow) 021 ■ 



where r = r{s,t) G R^, t is the time and s is arc length parameter along 
the vortex filament, finally a and f3 are assumed be constant (a > 0). The 
case (3 = corresponds to Localized Induction Equation of da Rios, closely 
associated with the cubic nonlinear Schrodinger equation and continuum 
Heisenberg ferromagnt model (see JI]). The case a = is associated with the 
modified Korteweg-de Vries equation. Usually a, (3 are treated as parameters 
of the theoretical model and they have to follow from experimental data. 

The exact formulas for multi-soliton solutions can be obtained using dif- 
ferent methods (compare jSl El CDl 1121 • follow the soliton surfaces ap- 
proach fn\ 11^1 as most convenient. The one-soliton solution is obtained 
without difficulties IT^ but two-soliton solution looks much more com- 
plicated (compare O CHI)- An essential simplification of the two-sohton 
formula is given in fUl El • 

However, the following statement written more than 16 years ago seems 
to be still accepted: "The shapes of the 'Hasimoto soliton' seem adequate to 
fit experimental data, but the ratio between group velocity and phase velocity 
for the wave is a constant (2.0) in LIA, whereas the experiments show a 
dependence on the ratio of wavelength to core size. Also unexplained is a 



r,t = a(r,sX 





(1) 
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large phase advance on collision of two experimental 'solitons'. An intriguing 
attempt to go beyond LIA as originally conceived was presented by Fukumoto 
& Miyazaki. Their proposal to include an arbitrary axial flow (■ ■ ■ ) leaves 
the shapes of the solitons unchanged, changes the velocity of the propagation, 
but, unfortunately, does nothing for the troublesome phase jump." [H]. 

The main goal of the present paper is to show that the large phase jumps 
are not longer "troublesome" and are in good agreement with the theory. 

In order to describe the interaction of two solitons on a single vortex 
filament we use the following exact formula for two-soliton solution: 



^ ( di (e^2 cosPi+ + e-^^ cosPf + rfa fe'^i C0SP2+ + e~^^ cosP^) \ 
di (e'^a sinPi+ + e''^^ sinPf j + ^2 (e'^' sinP^ + e"^^ sinPg") 
V rf+ sinlXQi + Q2) + rf- sinh(Qi - Q2) + sin(Pi - P2) - 2sK j 

(2) 

where 

K = coshQi cosh(52 coshA + sinhQi sinh(52 sinhA + sinhA cos(Pi — P2) , 

P± = Pfc±5, (/e = l,2), 

Qk = 2s ImAfc + 2t Imook + Qko , 

Pfc = —2s ReAfc — 2t Recj^ + P^o , 



d± = {di ± d2)e , do = di sin^i = — ^2 sin52 , d^ 



I A, 
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(by the way: changing the function co'(A) we get other equations of a large 
hierarchy of integrable systems, compare [18 ), A^ are complex parameters, 
and, finally, Qko, Pko, ^1,62, A are real constants. 

The constants Qko, Pko are arbitrary but not very important, they just 
fix the initial positions of the solitons. The constants 61,62, A depend on A^ 
as follows: 

gA+j<5i (Al - A2)A2 gA+jfe (A2 - Ai)Ai 

(Ai — A2)A2 (A2 — Ai)Ai 
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Two complex parameters Ai, A2 can be replaced by four real parameters 
01,02,61,62 (where = + ibk), or, which is even more convenient, by 
Cfc, dk expressed by: 

_ Qfc , _ bk _ Ck , _ dk . . 

One can also use the parameter Tk (compare |^) responsible for the shape 
of the soliton: Tk := Ck/dk = cik/bk- The parameter A depends on Ck,dk as 
follows: 

, ^1, (c^i - d2y + (ci - c^y , . 

2 (rfi + + (Ci - C2)2 

The two-soliton formula (j21) is obtained in the framework of the soliton sur- 
faces approach 16, 17j. Its form is relatively simple (especially when com- 
pared with other, equivalent, formulas, see jHl UHl Cll CS] ) and implies a lot 
of interesting consequences. Formulas for A^-soliton solution can be found 
for instance in [TH] (see also ^ ITIH ITTj). 

Let us first discuss the one-soliton solution (in the case ci;(A) = — 2A^ 
found by Hasimoto jS], extended for an arbitrary analytic ti;(A) by Sym [T^): 



ri 



/ 



V 



coshQi 



/ cosPi 

sinPi I . (6) 
y sinhQi 



Physical characteristics of the single soliton solution ri can be computed 
in the standard way. The amplitude of the Hasimoto vortex (i.e., its maximal 
distance from the z-axis) is given simply by di. The maximum of the wave 
envelope {Qi = 0) performs a helical movement and its z-coordinate (zmaxi) 
moves at a constant speed (group velocity vf): 

vf = W.i(t) = = -4(aai + PiSaj - bj)) . (7) 

ImAi 

The solitons can be classified into 3 classes: lump solitons (|ai| > \bi\), 
cusp solitons (|ai| = |6i|) and loop solitons (|ai| < \bi\) [Hlj. The names 
correspond to the shape of the curve z = z{p) (where p = \^x'^ + y'^). In 
experiments only the lump solitons were observed. If the vortex filament has 
some loops, then the localized induction approximation is not longer valid 
and these additional interactions probably can destroy the loops. However, 
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we remark that in the case of smoke rings (where such interactions are much 
more important) the locahzed induction equation was successfully apphed 
yielding some patterns which can be observed experimentally ([IDl, compare 

EH). 

The size of the lump soliton can be estimated by its half-width. The half- 
width Di is defined as a difference between the two values of z corresponding 
to Qi satisfying cosh Qi = 2: 

= di (^-V3 +(^1+ (^^) ln(2 + V3)^ (8) 

Another possibility is to estimate the wavelength of the soliton (i.e., the 
distance along the z-axis between neighboring wavecrests) 

Ai ^ - . 9 

ai 

Not entering into technical details we mention that more precise computa- 
tions lead to the result Ai = vr/ai + di(tanhQ[ — tanhQ'/), where Q'l and Q'l 
depend on s and t. The wavelength Ai is not constant along the filament and 
depends on t. In experiments usually the central wavelength is measured. 
The wavelength Ai equals vr/ai with an accuracy 2(ii/Ai which is less than 
l/vr. Usually, for di -C Ai, the accuracy is much better and, therefore, Ai 
is treated as a constant value given by 0- Note that always Ai > 2ndi. 
Obviously, the formulas ((Tj),®,© are valid for the second soliton {k = 2) as 
well. 

The formula Q enables us to compute the asymptotic behavior of r2. 
Assuming vf ^ ff we consider the limit Q2 ±cx3 (or, more precisely, 
I Qi I < I Q2 I). Thus 



Q2^±oo 

r2 — ' 




di 



cosh(Qi ± A) 



/ CO<Pi±5i) 

sin(Pi±5i) I (10) 
\ sinh(Qi ± A) 



The shape and the velocity of the soliton do not change during the interac- 
tion. The only result of the interaction is the phase shift. In fact we have 
two phase shifts: the shift A^^ along 63 axis and the shift of the angular 
variable Pi. The first one is much more important and can be measured in 
experiments. Analogical considerations are valid in the limit Qi ±00. 
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Let us consider the case of the head-on coUision of two sohtons (this 
case was observed in experiments jS]). We assume ff < 0, ff > (head- 
on colhsion) and also 6^ > (a technical assumption which can be made 
without loss of the generality). We are going to compute the phase shift of 
the first soliton (indexed by 1) after the interaction with the second soliton. 
If t ^ ±oo, then Q2 —>■ ±00 and the 2;-coordinate of the envelope maximum 
for the first soliton in asymptotic regions {Q2 ±00) can be easily derived 
from (Uni): 

(11) 

The phase shift Af^ := {z^^xi " -2^maxi)sgn(f f) and, taking into account 
vf < 0, we have: 

^f^.2,,-^^-2,,^i±^, J'j^f;if^-^% . (12) 

Oi 2di (di - d2)^ + (ci - C2)^ 

Analogical considerations lead to the similar formula for the phase shift of 
the second soliton: Ag'* = —2di — A/62- 

To compute the phase velocity of a soliton we have to determine positions 
of the individual wave peaks. Considering the projection of the soliton on the 
xz-plane we have to solve (at least implicitly) the equation ^ ( cosh o'fc ) ~ 
The resulting velocity of a wavecrest depends on its position Qi. 

,f = 2a'I^ + imi - 4) - HM(^±|M L _ n \ 

ak ctfe cosh Qk \ [ai + hi) cosh Qk j 

(13) 

In the analysis of solitons on a vortex filament the phase velocity has been 
always assumed to be constant sufficiently far from the envelope maximum 
[Qk 00). However in experiments usually the central region of the kink 
(i.e., Qk ~ 0) can be studied with the best accuracy. Anyway, the theoretical 
treatment of the phase velocity is relatively difficult and have to be made with 
care. Experimental results concerning the phase velocity are also father 
crude" |Zll8|. 

Another quantity estimated from experimental observation is the rate of 
rotation of the filament maximum (fi^), namely -Pfc(i)|Qfc=o = ^fc^ + const, 
where 

^]fc = -4(a + 4/5afc)(a^ + 6^) . (14) 
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To compare the presented theory with experimental results we recall 
that the two-soliton solution is parameterized by only 4 essential param- 
eters 01,61,02,62 (the parameters Qko,Pko, related to initial data, are not 
important for our analysis) and two more parameters {a, (3) characterize the 
physical model. In experiments one can measure more or less precisely much 
more quantities: amplitudes (rf^), torsions, velocities {vl,vl'^), widths {D^, 
Afc), angular velocities {flk), phase shifts (A^'^ and angular phase shifts) etc. 

We choose dk and Dk (or A^) as fundamental quantities. They can be 
measured with reasonable accuracy. They are very convenient in theoretical 
analysis beaucause from (jH)) we can express Ck in terms of dk and Dk. 



Ck = ±dk 



^ ln(2 + VS) 



(15) 



and all other characteristics of the solution are expressed by Ck, dk (by virtue 
of (jH) we can use the variables Ck, dk instead of a^, bk). From © we can also 
express Ck by A^: 

Ck Afc 



dk 2ndk \ \27idi. 



i-k' 

where > is assumed. Combining (fT3j) and (fT^ we get 



(16) 



Dk = (2Mik [^k + - Ij - l-73j dk , (17) 

where ^k '■= ^k/i'^'^dk) (we recall that ^k > !)• 

According to the typical experimental situation we assume di ~ ^2 and 
denote d := di. Then 

Af ^ (-2 + 0.66A)d + 0.38AL)fe (18) 

where A := ln(l + Ad^{ci - 02)"^). 

From dTTj) it follows Dk > A^ for Ak/dk > 10 and Dk > 2Ak for Ak/dk > 
17. Moreover, A > 1 if |ci — C2I < 1.53d, which can be estimated as |Ai — 
A2I < 5d. Then A^^ ~ O.ADk but this relation can easily "improved" (e.g., 
A^^ ~ Dk) if solitons are wider and more similar to each other. We arrive at 
the conclusion that head-on collision of lump solitons of similar size has to 
result in large phase shifts of both solitons. 
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The experimental observations of the sohton interactions have shown that 
"the forward advance (■ ■ ■ ) is of the order of the central wavelength of the 
isolated kink wave" (fH'], compare Fig. 4). The analysis presented above is in 
perfect agreement with this observation. For A 3 we have a very interesting 
coincidence: the phase shift is equal to the halfwidth of the soliton. 

We stress a very important point: the formulas (jl2p . (jl8j) are valid for any 
a, (3. Therefore our conclusion does not depend on the parameters of the 
physical model (provided that we confine ourselves to the localized induc- 
tion approximation of course), actually we can take even more general 
integrable equation described by an arbitrary function ci;(A). 

We propose the following procedure of comparing the presented theory 
with experimental data. First, taking and close to experimental values 
we use dSJ to compute ci, C2 (see (fT3jl ). In order to get a large phase shift we 
need |ci — C2I to be rather small, therefore ci and C2 should be of the same 
sign (this sign is related to the sign of a, which in turn depends on the the 
direction of the vorticity field along the filament, we take a > 8j). In any 
case we can use (fT2|) to choose right signs of Cfc. Then, we use ((Tj) to compute 
, ff • Thus we have a system of two linear algebraic equations for a and 
(3. The system is invariant with respect to the transformation: a — —a, 
Cfc — > — Cfc which means that an appropriate signs of Ck yield a > 0. Finally, 
having all six parameters, we may compare other quantities (e.g., rotation 
rates and phase velocities of the solitons) with experiment. Some quantities 
(e.g., phase shifts) need only 4 parameters to be computed (see above). 

Konno and Ichikawa JU]; trying to fit their two-soliton solution to the ex- 
perimental results of [H], explained very well amplitudes and group velocities 
of both solitons, to some extent their widths, and also the angular velocity 
and the phase velocity of one of the solitons (the corresponding experimental 
data for the second soliton were not reported in [H|). We are able to explain 
(with reasonable accuracy) the same data but also both phase shifts. 

Experimental data from [3 |H] (in units based on centimeter, second and 
radian) for the first soliton read as follows: di = 0.49, vf = —33, Ai = 6.8, 
Qi = —2.4, = —30, the torsion ti = 2ai = 0.89. The second soliton is 
described with lower precision: ^2 ~ (ii, ff ~ —vf, 2A2 ~ Ai. Finally, Af 
is "of the same order" as A^. My rough estimate, based on [HJ (Fig. 4), is: 
vi = 30, A2 ~ 3.5, Af ~ 5, Af ~ 1.5. 

Performing the above procedure we get the following parameters: ai = 
0.46, bi = 0.11, a2 = 0.78, 62 = 0.28. Then we assume vf = -33, v^ = 30 
and finally we get: di = 0.49, ^2 = 0.41, Ai = 6.8, A2 = 4.0, vf = -30.5, 
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Qi = -2.5, Af = 5.2, Af = 1.4. 

The agreement with experimental data is quite good (especially for the 
first soliton). Trying to get d2 closer to the value 0.49 we noticed an essential 
decrease of Af^ or increase of fli. Taking d2 = 0.41 is a kind of compromise. 

In this paper I revisited the Localized Induction Equation to show that, in 
spite of popular views, that model explains interactions of solitons observed 
in experiments including the existence of the large phase shifts. A little 
bit surprising is that my result is equally true for the Localized Induction 
Equation with and without axial flow. 

The theoretical investigation of the Localized Induction Equation, carried 
out in many directions (compare CHI 113 1221 12SI I21| ) , usually are not easy 
to be used by experimentahsts. I hope that the results of the present paper 
enable much better analysis of existing experimental results (which seem to 
be published only fragmentarily) and will prompt further experimental work 
in this field. Of special interest should be any observations of the interaction 
of vortex solitons in the superfiuid helium because such vortices are really 
thin which is assumed in the localized induction approximation. 
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